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■ We review the nucleon's twist-3 polarized structure functions from the viewpoint of gauge invariant, 
I I nonlocal light-cone operators in QCD. We discuss a systematic treatment of the polarized structure 

^ ^1 functions and the corresponding parton distribution functions. We emphasize unique features of higher 

QJ I twist distributions, and the role of the QCD equations of motion to derive their anomalous dimensions 

. for Q^-evolution. 

a^ : 

, I 1. Introduction 

>: 

, In the last ten years, great progress has been made both theoretically and experimentally in hadron 

' spin physics. Furthermore, in conjunction with new projects like the "RHIC spin project" , "polarized 

m; HERA", etc., we are now in a position to obtain more information on the spin structure of nucleons. 

■ In the case of spin-dependent processes, it has been known that the twist-3 contributions can be 
I measured as leading effects in certain asymmetries. Since the higher twist effects describe the coherent 

^\ • quark-gluon behavior inside the nucleon, we could obtain information on the correlation of quarks and 

I gluons beyond the (QCD) parton model (leading twist effect). 

^ • In this talk, we will concentrate on the twist-3 polarized structure functions ||l|. We present a systematic 

treatment of these functions in terms of gauge invariant, nonlocal light-cone operators, and summarize 
the recent theoretical progress on the QCD evolutions M. 



2. Classification of twist-3 structure functions 

To describe a variety of high-energy processes in a universal language, it is desirable to have a definition 



Oh! 
> 

X: 

JH , of parton distribution functions based on the operators in QCD. The traditional approach relies on the 

operator product expansion (OPE), but it can be applied only to a limited class of processes. This calls 
for an approach based on the factorization as a generalization of the OPE. 

2.1. Parton distributions 

It is instructive, at first, to consider the deep inelastic scattering (DIS) in free field theories. The 
result provides a guide to operator definitions of the structure functions (parton distributions). As is well 
known, the relevant quantity for the DIS is the hadronic tensor: 

W,. = J ^e-'^'^-{PS\[MO) , Mz)]\PS) . (1) 

Here \PS) is the hadron state with momentum P and spin S, (Q^ = —q^) is the momentum transferred 
to the hadron, and is the hadron's electromagnetic (weak) current. In the Bjorken limit, the hadronic 
tensor (|l]) is governed by the behavior of the current product near the light-cone (light-cone dominance). 
In the free field theory, the current product can be calculated to be 

[J^(0) , J,{z)] = ~ (a"A(z)) [5^a.a^^y(0, z) ~ le^c^uMliO, z)], (2) 
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up to irrelevant terms to the DIS. A{z) is the commutator function, 5^ 
and the nonlocal operators are defined by 

W^(0,z) = V;(0)7''VW - ^W7"^(0), (3) 
Ul(Q,z) = V5(0)7W'W + V^W7V^(0). 

The singularity of A(z) (or more generally the light cone dominance) selects the integration region 
with 77 ^ when we write as z^ = rjp'^ + Aw^ + z'^ using two auxiliary light-like vectors 

p2 = ^ , p-w ^ 1. (We take P'' = p^' + {^'P/2)w^' and = -xp^' + {Q^ /2x + M^x/2) w^" + 0{1/Q^) 
with the hadron mass M and the Bjorken variable x = Q'^ /{2P ■ q).) Therefore, we can expand the matrix 
element of (^ in powers of the deviation from the light-cone, — Aw^, and approximate it as 

(P5|Wf,^(0,z)|P5) ^ {PS\U^^AOAw)\PS). 

By parameterizing the first term of (^ as 

{PS\4>{0)r^{Xw)\PS) ^ 2p''q{X) + 2w''M^h{X), (4) 

and substituting it into (|^), the unpolarized structure function (quark distribution) is defined as the 
Fourier transform of q{X) into the Bjorken-x space, 

l(^) = / ^e^^^PSim jh^iXw)\PS). (5) 

The second term of (^ is the twist-4 contribution. Other terms of (||) can be treated similarly. These 
results yield the factorization into the c-number coefficients (short-distance part) and the parton distri- 
butions (long-distance part). 

In the presence of QCD interaction, the above results are modified in two respects. First, the higher 
order interactions produce the logarithmic (ln(Q^/^^)) corrections to A(z), and correspondingly the 
parton distribution functions acquire a dependence on the renormalization scale /i. Second, the coupling 
of the "longitudinal" gluons replaces A(2;) as Q 

A(z)^A(z)[0,z], 

where 

[y, z] = Pexp {^aj^d^iy - z)t.A''{ty + (1 - t)z) 

is the path-ordered link operator connecting the points z^ and along the straight line. By absorbing 
this factor, the nonlocal light-cone operators which define the parton distribution functions now preserve 
gauge invar iance. 

2.2. Twist-3 quark distributions 

Now, it is not difficult to define the parton distribution functions in QCD We consider the following 
quantity involving the nucleon matrix element: 



2^' 



-e 



iXx 



(F5|V'(0)[0,Aw]r?A(Aw)|F5) . (6) 



Here F is a generic Dirac matrix, and the link operator [0, Xw] makes the operators gauge invariant. 
Equation (^ defines the distribution function for a quark with momentum k ■ w = xP ■ w = x. F can be 
any Dirac matrix, depending on which hard process is considered. An important observation made in |^ 
is that one can generate all quark distribution functions up to twist-4 by substituting all the possible F. 
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By decomposing into independent tensor structures, one finds nine independent quark distribution 
functions associating with each tensor structure We hst here only the twist-3 polarized distribution 
functions. 

J ^e^^-{PS\m[0,Xw]r^rJj{Xw)\PS) = 2gT{x,^?)Sl + ■■■ , 
j ^e*^^(F5|V;(0)[0,Aw]a'"'i75V'(Au;)|F5) = 2 hL{x, ^?)M{p^'w'' - p''w^'){S ■ w) + ■ ■ ■ , 

where 5^ = S*!^ + S"^ with Sww^S-w. 

Three comments are in order (1) Since the definition of the twist here is based on the simple 
power counting with respect to l/Q, there is a slight mismatch with the conventional definition as 
"dimension minus spin" of the relevant operators. Actually both distributions qt and contain the 
contributions from the twist-2 pieces. (2) qt and Hl have different chiralities ("even" and "odd"). (3) 
The light-cone quantization formalism ^ is a conceptually useful approach to the twist counting. In 
this approach, quark fields are decomposed into "good" and "bad" components. The good component 
represents an independent degree of freedom; the bad components are not dynamically independent 
and can be reexpressed by a coherent quark-gluon pair. Therefore, only the twist-2 distributions are 
literally the "distributions" and correspond to the parton model, and the higher twist distributions are 
the multiparton (quark-gluon) correlations. 

2.3. Twist-3 gluon distributions 

Our analysis can be extended to the gluon distribution functions. The gauge-invariant definition of the 
gluon distribution functions is provided by 0] 

/rj\ J. 
— e^^^(F5|u;„G"^(0)[0, \w]w pG^'' {\w)\P S) = --Qsrix, y?)ie^''"^^ S^^wp + • • • , 

where G^'^ is the field strength tensor, and we show the twist-3 term explicitly. The gluon distributions 
mix through renormalization with the flavor singlet chiral-even quark distributions. On the other hand, 
there exists no gluon distributions that mix with the chiral-odd quark distributions. 

2.4. Twist-3 three-particle distributions 

Coherent many-particle contents of the nucleon are described by multiparton distribution functions. 
The twist-3 quark-gluon correlation functions are defined as 

/ ^^e'^'"'*'^^''~'^^'5|^(0)r[0,H5G^''(C^i')[Cw,Azi;]V(AuO|P5). (7) 

Similarly to the quark distributions, one can define the multiparton distributions by considering possible 
Dirac matrices for F. The treatment here can be extended to the case of three-gluon correlation functions, 
which are relevant to the (singlet) quark distribution grix) and the gluon distribution Gsrix). Although 
the extension is straightforward, we do not go into the details of the three-gluon nonlocal light-cone 
operators |H, 

3. QCD evolutions 

The calculation of evolutions for the higher twist terms is generally very complicated due to the 
presence of multiparton distributions. Although it is possible to generalize the DGLAP approach to 
the three-body case and calculate its evolution kernel, one is forced to use some particular techniques, 
e.g. the use of the light-like axial gauge etc. |l0| In this talk, we report an approach based on the 
"local composite operators" in a "covariant" gauge |ll| . The advantage of this approach is that all the 
relevant steps can be worked out based on the standard and familiar field theory techniques. Therefore, 
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the calculation is very straightforward and we can assess unambiguously the necessary information for 
the evolutions. 



3.1. Local operators 

By taking the moment of distribution functions, one can show that there is one-to-one correspondence 
between the nonlocal and the local operators. The simplest example is provided by the unpolarized 
twist- 2 quark distribution (|^): 

n-l ~j 

dxx^-^q{x) = t \ 5(A) I (PS'lVXO) ij\^,\w\il){\'w)\PS) 




= \ (^■5'|V'(0) ij {iw ■ Df'^ ^Q)\PS) . 

This equation implies that the moments of the parton distributions are given by the matrix elements of the 
corresponding gauge- invariant, local composite operators. For the twist-3 distributions, the situation is 
the same except that one should take into account the "double" moments for the three-body distributions. 

3.2. Operator relations 

The common feature for the higher twist operators is that there appear a set of local operators with 
the same quantum numbers, the number of which increases with spin (moment), and these operators are 
not all independent but related through the QCD equation of motion |l3|. 

One of the most convenient method to identify the relevant operators and a relation among them is 
to utilize the exact operator identities satisfied by the gauge invariant nonlocal operators. In the case of 
non-singlet part of gr, the following operator identity can be obtained by explicit differentiation 

d 



-^{Q)l"l^[Q,z]ij{z) - [a ^ ^l) 



dti^mQ. tz] / (^t ~ gG''P{tz)zp - ^gG^P{tz)zpj [tz, z]^{z) 



-f 2m,V;(0)75a'^''zp[0,z]V'(z) 

+ ^i0h5a^''zp[0,z]E^{z)-^{0)'E^^a^Pzp[0,z],p{z). 

Here E = i p — rriq is the equation of motion operator and represents the quark mass generically. 
This identity is exact through twist-3. Taylor expanding the above equation around z^ = 0, we easily 
identify the local operators, which contribute to the moment oi qt, and also the relation among them. 

Kl = :^{Vl-Vn-l-l+Ul+Un-i-l), 

Zn 

Zn 

where S symmetrizes over fii and A antisymmetrizes a and /ii, and 

Vi = -i'^gSipD''^ ■ ■ ■ G^^"' ■ ■ ■ D^'^-^J^'"-^J5'lP , 

Ui = -i"-^gS'4jD''' ■ ■ ■ G""' ■ ■ ■ D^—^j^^-'iP . 

The Rnj is related to the double moment of quark-gluon three-body correlation functions (j^) with 
r = 7^1, 7^1 75. These operators satisfy, 

n 1 "^^ 

n ^-^ 

1=1 
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For the singlet part, the situation becomes much more compHcated. We have additional operators 
that are made of two or three gluon fields, as well as the BRST invariant, gauge variant operators. The 
explicit expressions for the operators and the derivation of the identity satisfied by them are found in 
Ref. [|. 

For other distributions of twist-3 and higher, we can follow exactly the same procedure. In the case of 
e.g. ft-L, from the quantity, 

one can generate towers of identities for the corresponding local operators |[ Q . 
3.3. Calculations 

Since we identify a complete set of operators, the next step is to calculate the anomalous dimensions 
(renormalization constants) for these operators to derive the evolution. For this task, the following 
theorem (|l5| is very helpful. This theorem tells us that in non-abelian gauge theory, three kinds of 
operators will mix under the renormalization of gauge invariant operators. (I)the gauge invariant oper- 
ator itself. (Il)thc BRST invariant operators. (Ill)the operators which are proportional to the equation 
of motion (EOM). These EOM operators involve both BRST "invariant" and "variant" ones ||ll|, 
Although the physical matrix elements of the BRST invariant and EOM operators vanish Q , it is nec- 
essary to consider these operators to complete the renormalization, in other words to obtain the correct 
renormalization constants for the gauge invariant operators. 

At the lowest twist level, the above complicated operator mixing does not come into play, because there 
exists neither an EOM nor BRST invariant operator of twist-2. The EOM as well as the BRST invariant 
operators always have smaller spin by at least one unit than the possible highest spin operators of the 
same dimension. At the higher twist (> 3) level, both EOM and BRST invariant operators participate 
in the renormalization mixing with the gauge invariant operators. This is a characteristic feature of the 
renormalization for the higher twist operators. 

The explicit calculations at the one-loop level appear in Refs. for gr and in Ref. |l| for hj,. 



4. Conclusions 



We have surveyed the twist-3 polarized structure functions from a viewpoint of gauge invariant nonlocal 
operators. One of the key points of our approach is that we preserved maximal (gauge (BRST) and 
Lorentz) symmetries of the theory at every step of investigation. Another point is use of the one-to-one 
correspondence between the nonlocal and the local operator approaches to investigate the Q^-evolution. 
We have clarified the role of the QCD equations of motion to reveal the interrelation between the different 
twist-3 distributions and to derive their Q^-evolution. 

The polarized twist-3 distribution functions and hj^ are measurable as leading effects and are 
of special interest. However, the phenomenology of twist-3 distributions is not as straightforward as 
that of twist-2 distributions due to the fact that the twist-3 distributions are essentially three-particle 
correlations: We need to determine many parameters to make precise predictions. In this respect, it is 
pointed out that a drastic and universal simplification of all twist-3 distributions occurs in the Nc oo 
limit Q. It would be a good news to solve the above problem for all practical purposes. 

We hope that various kinds of new experiments and more detailed theoretical investigations will be 
able to clarify not only perturbative but also nonpcrturbative aspects of QCD related to hadron spin 
physics. 
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